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ABSTRACT 

We present predictions for the structural and photometric properties of early-type 
ON \ galaxies in the cold dark matter cosmology (ACDM) from the published semi- 

analytical galaxy formation models of Baugh et al. and Bower et al. These calculations 
were made with the GALFDRM code, which tracks the evolution of the disc and bulge 
components of a galaxy, using a self-consistent model to compute the scalelengths. 
The sizes of galactic discs are determined by the conservation of the angular momen- 
^j. ' turn of cooling gas. The sizes of merger remnants are computed by applying the virial 

ly-j , theorem and conserving the binding energy of the progenitors and their orbital en- 

OO ' ergy. There are a number of important differences between the two galaxy formation 

models. To suppress the overproduction of bright galaxies, the Bower et al. model 
^sO . employs AGN heating to stifle gas cooling, whereas the Baugh et al. model invokes 

a superwind which ejects cooled gas. Also, in the Baugh et al. model a top-heavy 
(— i \ stellar initial mass function is adopted in starbursts. We compare the model predic- 

O tions with observational results derived from the SDSS. The model enjoys a number of 

I 1 notable successes, such as giving reasonable reproductions of the local Faber- Jackson 

relation (velocity dispersion-luminosity), the velocity dispersion-age relation, and the 
fundamental plane relating the luminosity, velocity dispersion and effective radius of 
spheroids. These achievements are all the more remarkable when one bears in mind 
that none of the parameters have been adjusted to refine the model predictions. We 
study how the residuals around the fundamental plane relation depend on galaxy 
properties. We examine in detail the physical ingredients of the calculation of galaxy 
sizes in GALF0RM, showing which components have the most influence over our results. 
We also study the evolution of the scaling relations with redshift. However, in spite 
of the successes, there are some important disagreements between the predictions of 
the model and observations: the brightest model spheroids have effective radii smaller 
than observed and the zero-point of the fundamental plane shows little or no evolution 
with redshift in the model. 
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1 INTRODUCTION 

Remarkably tight correlations exist between the struc- 
tural and photometric properties of galaxies. Across the 
Hubble sequence there is a strong dependence of lu- 
minosity on either the rotation speed of galactic discs 



Tullv fc Fished 1977ft or the velocity dispersion of spheroids 
Faber fc Jackson! 1 19761 ) . Other scaling relations observed 



for early-ty pe galaxies include those betw een colour and 
magnitude (ISandage fc Visvanathar] ll978al fbT). colour and 
velocit y dispersion dBernardi et ~ 20051 ). radius and lumi- 
nosity (|Sandage fc Perelmuterl Il990ft ~ and radius and sur- 
face brightness (|Kormendvl 19771) . Some of these correlations 
can be combined into a "fundamental plane" which con- 
nects the effective radii, velocity dispersions and luminosi- 
ties of ellipticals (iDiorgovski fc Davis! 1 19871 ; iDressler et al.l 
ll987l : lBernardi et alj|2003al ). 



The existence of these scaling relations and their tight- 
ness encode clues about the formation and evolution of ellip- 
tical galaxies. For example, the existence of a fundamental 
plane can be understood by applying the virial theorem to a 
gravitationally bound stellar spheroid in dynamical equilib- 
rium, after making the assumption that ellipticals of differ- 
ent sizes have the same structure (homology) and a constant 
mass-to-light ratio. The deviation of the observed funda- 
mental plane from this prediction can therefore be driven 
by variations in the mass-to-light ratio across the early-type 
population or by a non-uniformity of the structure of ellip- 
ticals, referred to as struct ural non-homology , or a combi- 
nation of these two effects (|Ciotti et al.lll99q ; iBertin et al.l 
120021 ; iTruiillo et al.ll2u04al ). At first sight, the small scatter 
around the observed correlations would appear to pose a 
challenge to hierarchical galaxy formation models, since the 
variety of merger histories in the models would lead one to 
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expect to a corresponding scatter in the properties of early- 
type galaxies. 

In this paper, we use a semi-analytical approach to 
model the properties of elliptical galaxies. Such models pre- 
dict the star formation and merge r histories of galaxies (for 
a review of this class of model, see lBaughll2006h . In general, 
two channels are considered for the formation of spheroids: 
galaxy mergers or secular evolution of the disc. We will de- 
scribe the formation of discs and bulges in more detail in 
Section 2. The first attempts to track the disc and bulge of 
a galaxy separately simply recorded the mass and luminos- 
ity in each component (Baugh, Cole & Frenk 1996; Kauff- 
mann 1996). The models have now progressed to a state 
where detailed predictions can be produced for the struc- 
tural pro perties of galaxie s in addition to their stellar pop- 
ulations. ICole et all (|2000h introduced a model for the sizes 
of the disc and spheroid components of galaxies: the size 
of a galactic disc is calculated by assuming conservation of 
the angular momentum of the gas as it cools and collapses 
in the halo; the size of the spheroid is derived by applying 
conservation of binding and orbital energy, and by applying 
the virial theorem to the merging galaxies. The Cole et al. 
scheme also takes into account the gravitational force of the 
dark matter and the reaction of the dark matter halo to the 
presence of the baryons (see Section 2 for further details) . 

Cole et al. tested their model for the sizes of galactic 
discs against the observed distributio n of disc scale lengths 
estimated bv lDe Jong fc Lacevl (|200Ch . and verified that the 
predictions of their fiducial model were in excellent agree- 
ment with the observations. 

Cole et al. did not test their prescription for predicting 
the size of galactic spheroids. This is the focus of our pa- 
per. In related studies, Gonzalez et al. (2007, in prep.) test 
the predictions of galaxy scale lengths for disks and bulges 
at z = and Coenda et al. (2007, in prep.), look at the 
evol ution of galaxy sizes . 

lHatton et all (|2003l ) used a similar scheme to that out- 
lined by Cole et al. to compute the sizes of spheroids in 
the GALICS model. However, these authors adopted a less 
realistic model for the scale size of galactic discs. In com- 
mon with many semi-analytical models, they assumed that 
the scale size of a disc is related to the virial radius of the 
host dark matter halo by ro = A-R200/2, where A is the di- 
mensionless spin parameter for the dark matter halo, which 
quantifies its angular momentum, and ifeoo is the halo virial 
radius. This ignores the self-gravity of the baryons and the 
contraction they produce in the central regions of the dark 
matter halo. Several papers have considered the origin of 
the fundamental plane and the role of gas-rich and gas-poor 
mergers using numeri cal simulations, which follow the dark 
matter and baryons dKobavashil 120051; iDekel fc Coxl [2006; 



iRobertson et al. I l2006l ; iBovlan-Kolchin et al.l l2006 ) 

In Section 2, we summarize our model, explaining the 
ingredients which are particularly pertinent to the formation 
of galactic spheroids. We first compare our predictions to the 
sample of earl y-type galaxies drawn from the Sloan Digital 
Sky Survey by iBernardi et all (|2005l ): the selection criteria 
are described in Section 3 and the comparisons between our 
model predictions and the data are given in Section 4. In Sec- 
tion 5, we explore the sensitivity of our model predictions 
to various physical ingredients of the models. The evolution 
with redshift of the model predictions for the scaling rela- 



tions is presented in Section 6. Our conclusions are given in 
Section 7. 



2 THE GALAXY FORMATION MODEL 

A comprehensive overview of the GALF0RM model of galaxy 
formation and the p hilosophy behind s emi-analytical mod- 
ellin g can be foun d in lCole et al.l (|2000h (see also the review 
by iBaughl I2006T). Important ex tensi ons to the model ar e 
described in iBenson et all (|2002r i and iBenson et all l|2003h . 
In this paper , we focus on the pr edictions of the model in- 
troduced by iBaugh et alj (|2005l ). These authors put for- 
ward the first fully consistent hierarchical galaxy forma- 
tion model which was able to explain the observed num- 
ber counts of sub-mm sources and the luminosity function 
of Lyman-break galaxies, at the same time as reproducing 
the observed properties of the low redshift galaxy popula- 
tion. In some insta nces, we also show predictions from the 
model described bv lBower et al.l (|2006h . which includes feed- 
back processes associated with the accretion of material onto 
a supermassive bla ck hole, using the m odel of black hole 
growth explained in lMalbon et al. I <|2006h . Predictions from 
the Bower et al. model c an be downloaded over the internet 
fsee lLemson et al"1l2006l ). 

We now give a brief overview of the GALF0RM model, 
referring the reader to the references given in the previous 
paragraph for further details. We then outline some of the 
differences between the Baugh et al. and Bower et al. models. 
Finally, we recap some of the ingredients of the model which 
are particularly important for determining the masses and 
sizes of galactic spheroids and discs. 

The aim of the GALF0RM model is to make an ab initio 
calculation of the formation and evolution of the galaxy pop- 
ulation, set in the context of a cosmological model in which 
structures in the dark matter form hierarchically through 
gravitational instability. The main physical processes which 
we incorporate into the model are the following: (i) The hi- 
erarchical merging and collapse of dark matter haloes, (ii) 
The radiative cooling of shock heated gas. (iii) Quiescent 
star formation in discs, (iv) Feedback processes driven by 
supernovae and by the accretion of material on to supermas- 
sive b lack holes in the case of the model of iBower et al.l 
(2006). (v) The effect of a photoionizing background of ra- 
diation on the interg alactic medium and on galaxy forma- 
tion (see [Benson^et^l] [20o3) . (vi) The chemical enrichment 
of the gas and stars, (vii) The decay of the orbits of galactic 
satellites due to dynamical friction. This can lead to merg- 
ers between galaxies which can trigger bursts of star forma- 
tion and a change in galaxy morphology (see the next sub- 
section). The model generates a star formation history and 
a galaxy merger history for a representative population of 
galaxies at any epoch. Each galaxy is split into two compo- 
nents, a disc and a bulge. The formation of these components 
is discussed in the next subsection. 

For completeness, we now give a lis t of the primary 
differe nces between the lBaugh et alj (|2005l ) and lBower et al.l 
(2006) models, even though the Baugh et al. model is the 
main focus of the comparisons presented in the paper: 



(i) Dark matter halo merger trees. The Bower et al. model 
utilizes merger histories for dark matter haloes drawn from 
the Millennium Simulation of the hierarchical clustering of 
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dark matter in a ACDM universe l|Springel et al.ll2005l ). 
The simulation covers a volume of 0.125/i" 3 Gpc 3 . The mass 
resolution of the trees extracted from the simulation is 
1.72 x lO lo ft _1 Af . The Baugh et al. model is not set in 
the context of an N-body simulation. Instead, a representa- 
tive sample of galaxies is constructed by considering a grid of 
dark halo masses. For each mass on the grid, realizations of 
merger trees are generated using a Monte Carlo algorithm 
based on extended Press-Schechter theory (see I Cole et aU 
2000). The mass resolution used in the Monte Carlo trees is 
a factor of t hree better than th at of the trees drawn from the 
Millennium. iHellv et al.l (|2003T ) compared model predictions 
obtained using Monte Carlo or N-body merger trees and 
reached the conclusion that the results are very similar for 
bright galaxies, with the two prescriptions giving divergent 
answers for the luminosity function at faint luminosities; 
for the resolution of the Millennium, the predictions for the 
luminosity function are robust down to around three mag- 
nitudes fainter than the break in the luminosity function. 

(ii) The suppression of bright galaxy formation. Hierarchi- 
cal models tend to produce too many bright galaxies unless 
some physical mechanism is invoked to regulate the forma- 
tion of massive galaxies. The Baugh et al. and Bower et al. 
models do this in different ways. Bau gh et al. adopt a su - 
perwind driven by star formation (see lBenson et all 12003). 
In this case, the wind drives cold gas out of the galactic disc 
and out of the gravitational potential well of the dark halo. 
The effectiveness of the wind depends upon the depth of the 
potential well. Such winds have been observed in massive 
galaxies, with inferred mass ejection ra tes which are com- 
parable to the star formation rate (e.g. IPettini et al1l200ll ; 
IWilman et al.l l2005). In the Bower et al. model, the cooling 
of gas in quasi-static hot gas haloes is suppressed, effectively 
cutting off the "fuel supply" for star formation. These are 
haloes in which the cooling time of the gas exceeds the free 
fall time within the halo. The cooling is quenched by the 
energy injected into the hot halo by the accretion of mass 
onto the central supermassive black hole in the galaxy. The 
growth of the black hol e is fo llowed using the model de- 
scribed bv lMalbon et all (120061 ). 

(iii) Quiescent star formation in discs. The scaling with 
redshift of the timescale for quiescent star formation in 
galactic discs is different in the two models. Both models 
allow the star formation timescale to depend upon some 
power of the circular velocity of the disc, and multiply this 
by an efficiency factor (eqn. 4.14 of Cole et al.). Baugh et al. 
assume that the efficiency factor is independent of redshift, 
whereas in the Bower et al. model this factor scales with the 
dynamical time of the galaxy. Dynamical times are shorter 
at high redshifts, so the quiescent star formation timescales 
are shorter in the Bower et al. model at high redshift than 
they are in the Baugh et al. model. 

(iv) Star formation in bursts. Globally, the two mod- 
els display somewhat different star formation densities in 
bursts. Furthermore bursts are triggered in different ways in 
the models. Baugh et al. only consider bursts resulting from 
galaxy mergers. A burst may accompany a major merger in 
which a galaxy accretes a satellite of a comparable mass to 
its own, or a minor merger in which a gas rich disc is hit by a 
much smaller satellite. Bower et al. consider this mode of ini- 
tiating star bursts, but also incorporate bursts which occur 
when discs become dynamically unstable. The need for this 



additional channel for bursts is driven by the need to build 
up black hole mass, so that cooling flows can be suppressed 
in massive haloes. One final difference to note between star 
bursts in the two models is that Baugh et al. invoke a flat 
initial mass function (IMF) for stars produced in a burst, 
whereas Bower et al. adopt a standard IMF (Bower et al. 
adopt a Kennicutt (1998) IMF in all modes of star forma- 
tion; Baugh et al. use the Kennicutt IMF in quiescent star 
formation) . 

(v) Background cosmology. Baugh et al. adopt the pa- 
rameters of the concordance ACDM model: matter density, 
Qo = 0.3, cosmological constant, Ao = 0.7, baryon density, 
ttb = 0.04 and a normalization of density fluctuations given 
by as = 0.93. Bower et al. use the cosmological parameters 
of the Millennium simulation, which are in better agreement 
with the constraints from the anisotropies in the cosmic mir- 
cowave background and large scale galaxy clustering (e.g. 
ISanchez et al]|2006l ): fi = 0.25, A = 0.75, Q b = 0.045 and 
as = 0.9. The power spectrum of density fluctuation used in 
the Millennium has somewhat more large scale power and 
less small scale power than that in the concordance ACDM 
model. 

One difference between the predictions of the two mod- 
els is the amount of star formation which takes place in star- 
bursts. Baugh et al. calculate that 30% of all star formation 
in their model takes place in starbursts driven by galaxy 
mergers. However, due to the high recycled fraction in star- 
bursts as a consequence of a top-heavy IMF, only 7% of this 
mass is locked up in long lived stars, and yet bulges account 
for around 50% of the global stellar mass in this model. The 
re-assembly of stars which were produced with the standard 
IMF used in quiescent star formation is therefore the pri- 
mary source of stellar mass in bulges (Baugh, Cole & Frenk 
1996; de Lucia et al. 2006; de Lucia & Blaizot 2006). In the 
Bower et al. model, the amount of star formation triggered 
by galaxy mergers is lower than in the Baugh et al. model 
for a combination of reasons: 1) Bower et al. use a quiescent 
star formation timescale which depends on the dynamical 
time, and is thus shorter at high redshift. Galactic disks at 
high redshift are therefore gas poor in the Bower et al. model 
compared with those in the Baugh etal model, so there is 
less fuel for the starbursts. 2) Minor mergers do not trigger 
bursts of star formation in the Bower et al. model, as is the 
case in the Baugh et al. model. In the Bower et al. model, 
disk instabilities account for over half of the star formation 
in bursts. Due to the choice of a standard IMF in all modes 
of star formation in this model, bursts are responsible for 
producing around a quarter of the mass in bulge. Again, 
the re-assembly of stellar mass made in galactic disks is the 
main source of spheroid stars. 

2.1 The formation of spheroids and discs 

We assume that gas initially settles into a rotationally sup- 
ported disc when it cools from the hot halo. This gas even- 
tually turns into stars in the quiescent star formation mode. 
The effective timescale on which the star formation takes 
place does not d epend upon t h e dyn amical time of the disc 
in the model of iBaugh et all (|2005h , but does have some 
dependence on the circular velocity of the disc. 

The formation of galactic spheroids takes place through 
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two channels: galaxy mergers and the instability of galactic 
discs. Baugh et al. only consider the galaxy merger mode of 
spheroid formation; Bower et al. consider both mechanisms. 

The consequences of a galaxy merger are characterized 
by the ratio, R = m sa t/m cen trai, of the mass of the accreted 
satellite galaxy (m sa t) to the mass of the primary or cen- 
tral galaxy in the halo (m cen trai), onto which the satellite 
is accreted. The mass ratio R is compared to two thresh- 
olds, /eiiip and /bu rst , to establish the sev erity of the merger 
(see, for example. iBournaud et al1l2005h . These thresholds 
are model parameters. If R > / c ni p = 0.3, then the galaxy 
merger is termed "major" . In the case of major mergers, the 
disc of the primary galaxy is destroyed. All stars are trans- 
ferred to the spheroid component and any cold gas present 
participates in a burst of star formation which adds stars to 
the spheroid. If the mass ratio of satellite to primary falls 
between the thresholds, i.e. /burst < R < f e m P , then the 
stellar disc of the primary survives and the stars from the 
accreted satellite are added to the spheroid. In this case, 
if the primary is also gas rich, that is if cold gas accounts 
for at least 75% of the total mass in the disc, then we as- 
sume that the accretion of the satellite induces an instability 
which drains the primary disc of cold gas, leading to a burst 
of star formation in the spheroid. 

In our model, a galax y can move in eith er direction 
along the Hubble sequence ijBaueh et al . 1996). The accre- 
tion of gas from the hot halo and subsequent quiescent star 
formation leads to a late-type (disc dominated) galaxy. A 
major merger between two such galaxies produces a descen- 
dent galaxy which jumps to the opposite end of the Hubble 
sequence, becoming an early- type galaxy (bulge dominated). 
Further accretion of cooling gas allows the galaxy to grow a 
new disc around its bulge, moving the galaxy back towards 
the late-type part of the sequence. 

In addition to the merger mode of spheroid production, 
Bower et al. also consider the secular production of bulges 
from discs which are unstable due to their strong self-gravity. 
This mode is most important in less massive galaxies. 



2.2 The scale lengths of the disc and bulge 
components of galaxies 

We assume that discs have an exponential profile, with a 
half-mass radius given by rdi sc , and bulges have a r 1 / 4 profile 
in projection, with a half mass in 3D given by rbuige- 

The scalelength of the disc is determined by the an- 
gular momentum of the halo gas, which arises due to the 
tidal torques which act during the formation of the halo. 
The angular momentum of the halo gas is quantified by 
the dimensionless spin parameter, A; this quantity is as- 
sumed to follow a log-normal dist ribution matchin g the re- 
sults of N-body simulations (see ICole et al.l 120001 . for de- 
tails) . We assume that the angular momentum of the gas is 
cons erved as it cools to f orm a rotationally supported disc 
(see lOkamoto et al. 2005, for a discussion of this assump- 
tion). 

Spheroids are formed in galaxy mergers or through disc 
instabilities as outlined in the previous subsection. The size 
of the spheroid resulting from a galaxy merger, r m , is deter- 
mined by applying the conser vation of energy and the virial 
theorem (see Section 4.4.2 of lCole et aHl2000h : 



(Mi + M 2 ) 



Mj 



T2 



/orbit MlM 2 
C n + T2 ' 



(1) 



where Mi represents the total mass (stellar, cold gas and 
dark matter) of one of the merging objects, within n , and 
the form factor c and the parameter /orbit are related to the 
self-binding energy and orbital energy by 



-Ebind 
-Eorbit 



GMf 



/orbit GM 1 M 2 

2 n + r 2 ' 



(2) 
(3) 



For simplicity, we adopt c = 0.5 and /orbit = 1. Later on, 
we explore the impact on our predictions of varying / or bit- 
Similar arguments are applied to calculate the scale size of 
the sph eroid which resul ts from an unstable disc (see Section 
4.4.3 of lCole et al.1l2000l) . 

Once the scale lengths of the disc and bulge compo- 
nents have been calculated as outlined above, we next take 
into account the selfgravity of the baryons and the contrac- 
tion of the dark matter halo in response to the gravity of 
the condensed baryons. New radii are computed for the disc 
and bulge by applying an adiabatic contraction of the disc , 
bulge and dark ma tter components (|Blumenthal et al1ll986l ; 
IJesseit et aTll2002l ). In the case of the disc, the total specific 
angular momentum is conserved. The bulge and dark matter 
halo are not rotationally supported. Nevertheless, it is use- 
ful to define an equivalent circular velocity using the veloc- 
ity dispersion of each of these components, and, using this, 
to define a quantity which we refer to as a pseudo-angular 
momentum. For the bulge the pseudo-angular momentum 
is given by: j b ui gc = r- b uig C V c (r b uig C ). This quantity is con- 
served during the adiabatic contraction. A similar quantity 
is conserved for the dark matter. 

In the case of the secular growth of spheroids, we again 
apply the conditions of virial equilibrium and energy conser- 
vation as in Eq. [1] defining the component 1 as the galactic 
bulge, Mi = Mbuig C , ri = rb u i gc , and the component 2 as the 
unstable disc, M2 = Mdi sc , ?"2 = ?"disc- After the calculation 
of the radius of the new bulge, we readjust adiabatically the 
spheroid and halo terms to reach the new equilibrium. 



3 THE SELECTION OF AN EARLY- TYPE 
SAMPLE 

We first compare our model predictions against the scaling 
relations and statistics of the sample of early-type galaxie s 
constructed from the SDSS bv lBernardi et all (|2003al . [2005l ). 
These authors measured relations between luminosity and 
various properties of early-type galaxies such as velocity dis- 
persion, effective r adius, effective mass, effectiv e density and 
surface brightness (|Bernardi et al.l i2003b, 20Q5|). The sample 
was also used to measure the lumin osity function of early- 
type galaxies ([Bernardi et al.ll2003bT ) , the fundamental plane 
(Bernardi et al l l2003d )~ and the colour-magnitude/colou r- 
velocity dispersion relations (Ber nardi et alj[2003dl 12005). 

For a complete description of the construction of the 
early-type sample from the SDSS, we refer the reader to 
the above papers. Below we give a summ ary of the selection 
criteria applied bv lBernardi et all (|2005l ) (hereafter Ber05). 
Galaxies are included in the sample if they have: 
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(i) Redshift z ^ 0.3, with a median z me d = 0.13. 

(ii) Apparent r-band Petrosian magnitudes in the range 
14.5 < r p 

ctro < 17.75. 

(iii) eclass < 0. The eclass value is a classification of 
the spectral type of a galaxy derived fro m a principal com- 
pone nt decomposition of its spectrum l|Connollv fc Szalavl 
1 19991 ). Ber05 chose negative values of eclass as this corre- 
sponds to spectra in which absorption lines dominate, char- 
acteristic of early-type galaxies. 

(iv) fracDev> 0.8, computed using the r— band image. 
The value of fracDev is an indicator of morphology. It is 
calculated in two steps, in turn. First, the best fit de Vau- 
couleurs and exponential profiles to the galaxy image are 
found. Second, using the scale lengths of the best fit pro- 
files found in step one, the best fit linear combination of the 
disc and bulge profiles is derived. The contribution of the de 
Vaucouleurs profile to this linear combination is the value 
of fracDev. 



We attempt to reproduce these selection criteria by im- 
posing the following conditions on GALFORM galaxies: 



(i) We generate a population of galaxies at an output red- 
shift of z — 0.13, the median redshift of the Ber05 sample. 
In our comparisons, we consider only Ber05 galaxies that lie 
within the redshift interval 0.11 < z ^ 0.15, close to the me- 
dian redshift. This additional selection in redshift reduces 
the size of the observational sample by a factor of ~ 4 to 
~ 11000 objects. 

(ii) We use total magnitudes to select a sample of model 
galaxies. We apply the same apparent magnitude limits 
which are used for the data. This is a reasonable approach 
as the difference between total an d Petrosian magnitu des is 
typically smaller than 0.2 mag <|Graham et al.ll2005h . We 
have also computed Petrosian magnitudes for our model 
galaxies and find that using the Petrosian magnitudes in 
place of total magnitudes does not make a significant differ- 
ence to our results. 

(iii) At present, GALFORM does not produce spectra with 
absorption line features. Therefore we cannot directly cal- 
culate a value for the spectral parameter eclass. Instead, 
we use the g — r colour which is more readily predicted for 
model galaxies. In Fig. [1] we use the SDSS DR4 to show 
that there is a good correlation between g — r colour and 
eclass. We retain galaxies with g — r > 0.8; Fig. Q] shows 
that more than 95% of the galaxies with a negative value 
for eclass are selected by this colour cut. 

(iv) We compute the value of fracdev in the a similar 
way as was done for the SDSS galaxies. We assume that the 
model bulges follow a de Vaucouleurs profile and the discs an 
exponential profile; these profiles describe the distribution 
of stellar mass, and so are independent of the passband. We 
adopt a cut on fracdev > 0.8. In Section 4.1, we explore the 
impact on our predictions of replacing the cut in the value of 
fracdev with a simple cut on the bulge-to-total luminosity 
ratio of the model galaxies. 

(v) Due to limitations of the SDSS data we also set a 
surface brightness threshold, \i e < 24.5 mag arcsec -2 . 
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Figure 1. The relation between g — r colour and eclass, for a 
sample of galaxies selected from SDSS DR4 in the redshift range 
0.12 < z < 0.15 and with Petrosian magnitudes in the interval 
14.5 < r pe tro < 17.75. The dotted line represen ts the selection 
applie d to the observational data, eclass < 0, bv lBernardi et all 
(2005): more than 95% of these galaxies have g — r > 0.8. 



4 RESULTS 

In this sec t ion, we compare the predictions of the 
Baug h~et al.l (|2005l ) model with observational data for early- 
type galaxies derived from the SDSS sample of Ber05. 



4.1 The luminosity function of early-type galaxies 

The luminosity function is perhaps the most fundamental 
statistical description of the galaxy population. Later on, we 
perform fits to the fundamental plane of early-type galax- 
ies in our model. The results for this fit are sensitive to the 
abundance of galaxies as a function of luminosity, so it is im- 
perative that the model reproduces the observed luminosity 
function closely. 

In Fig. [2] we compare the predictions of the GALFORM 
model for the r-band luminosity function of early-type 
galaxies with the estimate from the SDSS sample of 
Ber05. The luminosity function of SDSS early-types 
is well described by a Gaussian form: <fi(M) dM = 

i) 2 /{2o 2 )}, where 



M* + 



(j)*/\ / 2na 2 exp{-(M 

(0*,M*,(T,<2) = (1.99 x 10~ 3 Mpc- 3 , -21.15, 0.841, 0.85) 
respectively (note, Ber05 assume h = 0.7). The model 
predictions are in reasonably good agreement with the 
luminosity function estimated from the data. 

We also show, in Fig. [2] the effect of changing the cri- 
teria used to select early-type galaxies in the model. In 
our standard selection, the primary indicator of morphol- 
ogy is fracdev (see Section 3). We have also explored using 
the bulge-to-total luminosity ratio i n the r-band, B/T r , in 
place of fracdev l|Baugh et alj|l996h . The results for cuts of 
B/T r > 0.8 and B/T r > 0.5 are shown by the dotted and 
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Figure 2. The luminosity function of early-type galaxies at 
2 = 0.13. The solid black l ine shows a fit to the luminosity func- 
tion of the SDSS sample of Bernardi et al The results for 
the GALF0RM sample are plotted using grey lines. The solid grey 
line shows our standard early-type galaxy selection, as outlined 
in Section 3. The errorbars show Poisson errors due to the finite 
number of galaxies simulated. The other lines show how the lu- 
minosity function varies when, instead of using fracdev > 0.8, 
the bulge-to-total r-band luminosity ratio is used; the dashed 
line shows the results for B/T r > 0.5 and the dotted line for 
B/T r > 0.8. 



dashed lines respectively in Fig. [2] The luminosity function 
derived using B/TY > 0.8 is remarkably similar to the one 
obtained using fracdev > 0.8 (shown by the solid grey line). 

In our subsequent comparisons with the SDSS sample 
of early-types, we assign each model galaxy a weight which 
depends upon luminosity, such that the luminosity function 
of early-types in the model is forced to reproduce exactly 
the luminosity function of the data. 



4.2 The Faber-Jackson and cr-age relations 

The Faber-Jackson (hereafter FJ) relation was one of the 
first scaling relations to be discovered for early- type galax- 
ies (|Faber fc Jacksonll 19761 ). This relation indicates that lu- 
minosity is a strong function of velocity dispersion, a, with 
brighter early-types displaying larger velocity dispersions. 
Observational studies sugge st that this rela t ion i s approxi- 
mately given by Lex g 4 : lForbes fc Ponmanl i|l999l ). using a 
local sample of early-type galaxi es, found Lb oc a 3 ' 9 in the 
B-band, while lPahre et all (|l998h reported Lk oc a 4 " 1 in the 
K-band. These results also indicate that the FJ relation is 
essentially indepen dent of wavelength. In the case of SDSS 
early-type galaxies, iBernardi et alj (|2003bl ) confirmed these 
earlier results, finding L r oc a 3 ' 91 in the r-band at z = 0, 
with no significant differences in slope apparent in the g, i 
or z-bands. 

Fig. [3] shows the velocity dispersion-magnitude relation 



Figure 3. The Faber-Jackson relation between luminosity and 
velocity dispersion. The GALF0RM prediction is shown in grey and 
the hatched shaded region shows the relation for the Ber05 sam- 
ple. The shaded and hatched regions connect the 10 and 90 per- 
centiles. The one-dimensional velocity dispersion was calculated 
using am = (l.l/x/JOK^bulgei as explained in the text. 



predicted by GALF0RM. The one-dimensional velocity dis- 
persion is calculated using aio = (l-l/v / 3)lc j buigc, where 
bulge is the effective circul ar velocity of the bulge, which 
is as sumed to be isotropic l|Frenk et alj 1 19881 ; I Cole et ah! 
1 1994 ). The factor of 1.1 is an empirical correction which 
Cole et al. employed to map data for elliptical galaxies onto 
the Tully-Fisher relation for spiral galaxies. Fig. [3] shows 
that retaining this factor gives good agreement with the ob- 
served FJ relation. We find no change in the predictions if 
we consider, instead, the effective circular velocity of the 
disc and bulge combined. This is a consequence of our selec- 
tion which ensures that the model galaxies we consider are 
bulge dominated, as shown by Fig. [2] We find reasonably 
good agreement between model predictions and the FJ rela- 
tion observed for the Ber05 sample, albeit with a shallower 



slope, L r oc a 



(note, we plot a on the y-axis). Whilst 



the slope of the predicted FJ relation is formally at odds 
with that measured by Ber05, it is clear from Fig. 3 that 
the velocity dispersion of the model galaxies would change 
by relatively little even in the case of perfect agreement be- 
tween the observed and predicted slopes, given the limited 
range of magnitudes plotted. We find little dependence of 
the slope of the FJ relation on passband, in agreement with 
observations. 

The evolution of the FJ relation with redshift is plot- 
ted in Fig. [4] Here, we have chosen output redshifts in 
GALF0RM to match the median redshifts of the Ber05 red- 
shift subsamples: for SDSS galaxies with z < 0.08 we use 
Zmcd = 0.06 for the model galaxies, and for SDSS galaxies 
with z > 0.18 we set z mc( j = 0.20. Little evolution is ob- 
served in the zero-point of the FJ relation with redshift, a 
trend which is reproduced by the model predictions. The 
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Figure 4. The evolution of the Faber- Jackson relation for early- 
type galaxies. The upper panel shows the model predictions at 
redshifts z = 0.06, z = 0.13 and z = 0.20, which are the me- 
dian redshifts of the observational subsamples shown in the lower 
panel. Shaded regions connect the 10 and 90 percentiles of the 
distributions. 



shift to brighter magnitudes with increasing redshift is sim- 
ply a reflection of the fixed apparent magnitude limit of 
the SDSS. As we shall see in next section, this absence of 
evolution is actually expected at these redshifts due to the 
cancellation of two different evolutionary effects. 

Finally, in Fig. [5] we plot the luminosity- weighted age of 
the stellar population, computed in the r-band, as a function 
of the velocity dispersion. Some authors have argued that a 
correlation exists between these quantitie s, which has impli- 
cations for the scatter in the FJ relation (|Forbes fc Ponmar] 



Figure 5. The r-band luminosity-weighted age of the stellar 
population as a function of the velocity dispersion predicted by 
GALFORM. The model galaxies are at redshift z = 0.13 and are 
selected in a similar way to the observational sample of Ber05. 
The contours enclose 35%, 72%, 89% and 99% of the number 
density respectively, moving from black to grey. 



1 19991 ; iNelan et al]|2005l ). Fig. [5] reveals that velocity dis- 
persion and luminosity-weighted age are indeed correlated 
in the model, with galaxies which have larger velocity dis- 
persions also displaying older stellar populations. A linear 
fit to the model predictions shows that Age oc o -°- 58 ± 02 5 
wh ich is in exc e llent agreement with recent determination 
by INelan et~aH <|2005l ). who found Age oc cr - 59 ± 013 . Fur- 
thermore, we verify that the inc l usion of AGN feedback, 
as implemented by iBower et alj (120061 ). does not change 
this relation substantially, giving Age oc a - 51±0 - 03 ^ At first 
sight, these predictions seem to contradict those presented, 
for a different observational selection, by Nagashima et al. 
(2005b). However, it is important to note that, at least in 
the case of the model galaxies, the slope and scatter of the 
Age — a relation are very sensitive to the selection criteria 
applied. 



4.3 Radius-Luminosity Relation 

Another component of the fundamental plane of early-type 
galaxies is the relation between radius and luminosity; galax- 
ies with larger radii are more luminous. This was originally 
of interest fo r use in distance scale measure ments and cosmo- 
logical tests (|Sandage fc Perelmuter|[l990l ). Different studies 
indicate that this relation v aries slightly with wavelength. 
For example, S chade et al. I (Il997t) determined Lb on 



in the B-band and lPahre et al.l (U998T ) foun d Lk oc 



the K- band. For SDSS early-type galaxies, iBernardi et alj 
(|2003bl ^ reported L g oc r^ 50 in the g-band and L oc r^ 58 in 
the r, i and z-bands, which is consistent with the variation 
of the slope of this relation with wavelength suggested by 
the earlier determinations. 
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Figure 6. The relation between half-light radius and r-band 
magnitude for early- type galaxies. Again, the predictions for 
GALFORM galaxies are shown in gray and the hatched shaded dis- 
tribution represents the SDSS sample; in both cases the shading 
shows the 10 to 90 percentile range. 

We compare the predicted radius-luminosity relation 
with the Ber05 data in Fig. [6] The effective radius plot- 
ted here, r c , is the projected bulge half-light radius of the 
de Vaucouleurs law, which is related to the half-mass radius 
in 3D, r b , by: r c — r b /1.35. The model predictions do not 
change if we compute a composite half-mass (half-light) ra- 
dius from the disc and bulge components of the galaxy, as 
the model galaxies we consider are bulge dominated. The 
slope of the predicted radius-luminosity relation is flatter 
than is observed. The agreement between the model predic- 
tions and the observations is best at fainter magnitudes; the 
brightest early-type galaxies are predicted to be around a 
factor of three smaller in effective radius than is observed. 
Fig. [5] and Fig. [3] suggest that in the model, the high lumi- 
nosity early-type galaxies have a pseudo- angular momentum 
which is lower than would be inferred from the data (see 
§ I2.2l for the definition of the pseudo-angular momentum of 
a spheroid) . 

The evolution of the radius-luminosity relation with 
redshift is plotted in Fig. [7] We find no clear change in 
the slope of the radius- luminosity relation over this red- 
shift int erval, in agreement with the results of lBernardi et all 
(2003b|). We shall return to this point in section \§\ 

4.4 Effective Mass 

We can define an effective dynamical mass, Md yn , by set- 
ting A/dyn = r e a 2 /G. This differs from the true mass, M, 
because the definition of Md yn ignores any rotational sup- 
port and the flattening of galaxies. The difference between 
the two masses can be quantified by a correction term, £: 
M — £A4dyn- For a galaxy with T-type E0, there is no flat- 
tening or rotational support and so £ = 1. tn contrast, for 
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Figure 7. The evolution with redshift of the radius - r-band 
luminosity relation. The upper panel shows GALFORM galaxies at 
the median redshifts of the observational subsamples: z = 0.06, 
2 = 0.13 and z = 0.20 . The lower panel shows the results for the 
iBernardi et al ] j2005l) sample divided into volume-limited bins, 
as indicated by the legend. 

the case of an E6 galaxy, the true mass is al most twice as 
large as the dynamical mass, with £ ~ 1.9 fsee lBender et all 
1992, for details). 

In Fig.[8]we compare our prediction for the relation be- 
tween dyn amical mass and lumin osity with the observed re- 
sult for the lBernardi et alj 1)20051 ) sample. The figure reveals 
reasonable agreement between the model and the observa- 
tions for fainter galaxies. Brighter galaxies, in the model, 
have a somewhat lower dynamical mass than observed. As 
we noted when discussing the radius-luminosity relation, 
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Figure 8. The relation between dynamical mass and luminosity. 
The GALFORM data is represe nted in gray and the dark hatched 
shaded region represents the iBernardi et al.l l|2005h sample. The 
shading connects the 10 and 90 percentile values. 



high luminosity galaxies in the model display a lower specific 
pseudo-angular momentum than is observed, which trans- 
lates into a smaller dynamical mass. 



4.5 Fundamental Plane 

Observational studies indicate that early-type galaxies show 
tight correl ations between their kinematic and photometri c 
properties (|Diorgovski fc Davis! [l987l ; iDressler et all 1 19871 ). 
The remarkably small scatter about the so-called fundamen- 
tal plane connecting the effective radius, velocity dispersion 
and surface brightness of early types encodes information 
about the formation and evolution of these galaxies. 

The existence of a fundamental plane is expected if 
a stellar system obeys the virial theorem, which connects 
the kinetic and potential energies. The assumption of virial 
equilibrium gives a relation between the three-dimensional 
velocity dispersion, <T3d and the "gravitational" radius, r g , 
assuming that the system is gravitationally bound: 



2 

0"3D 



GM 



(4) 



This equation can be rewritten in terms of the central one 
dimensional velocity dispersion, ctid, and the effective radius 
of the galaxy, r c , 

GM 



where we have defined structural constants such that 



(5) 



based on the assumption that the population is homologous. 
The mean surface brightness within half-mass radius 



of a galaxy is I c = L/2nrl, where L is the total luminos- 
ity of the galaxy and the mean surface density is given by 
E c = M/2nr^. The ratio of the surface brightness to the 
surface density is equal to the mean mass-to-light ratio of 
the galaxy, within r c : E c // C = M/L. Using these definitions, 
Eq. [S] can be rearranged to give an expression for the fun- 
damental plane, 



r c = 



_2 
<?1D 



lo! 



2nG I e (M/L) ' 
;r c = 21oga + 0.4/ie +log(^r^v) 



log(M/L)+ 7 ,(6) 



where 7 is a constant whose value depends upon G and the 
choice of units. 

The observed plane is slightly different from the form 
predicted in Eq. [6] which follows by applying the virial 
theorem to a purely stellar galaxy without any dark mat- 
ter and assuming an h omologous population. For example, 
|j0rgensen et all (I1996I 1 ) found log r e = 1.24 log a + 0.328/i c + 
V. while IBernardi et ail (|2003cr ) obtained logr e = (1.49 ± 
0.05) logo- + (0.30 ± 0.01)/x e - (8.78 ± 0.02). 

The discrepancy between the theoretical prediction out- 
lined above and the obser vational result s is kno wn as the tilt 
of the fundamental plane. [Truiillo et al.l (|2004al ) argued that 
this tilt is due to a combination of effects: structural non- 
homology, which means a change in the surface brightness 
profile of the early-types with luminosity, and a variation in 
the mass-to-light ratio of the stellar populations with galaxy 
luminosity. 

The intrinsic thickness or scatter in the fundamental 
plane pose s another cha l lenge , and its interpretation is far 
from clear. [Forbes et all (1998) showed that the scatter was 
mainly due t o the age of the stellar population. However, 
iPahre et all (|l999h demonstrated that the position of the 
galaxy relative to the FP could not be entirely due to age 
or metallicity effects. 

To determine the location of the fundamental plane, we 
consider an orthogonal fit to the plane given by: 



log 



a log a ± b fi c + c, 



and determine the values of the coefficients a, b and c by 
minimizing the quantity 



S = 



N n 



a log a.i —bfi e 



i=l 



1 + a 2 



b 2 



Following this procedure, we obtain a fundamental plane 
for GALFORM given by log r c = (1.94 ± 0.01) logo" + (0.19 ± 
0.01)ai c - (7.54 ± 0.03) in the r-band. 

In Fig. [9] we plot the fundamental plane derived 
from GALFORM model galaxies, along with the data from 
IBernardi et"aH (|2005h in the same projection of the plane. 
Fig. [9] reveals reasonable agreement between the funda- 
mental plane predicted by GALFORM and the observational 
data: we can reproduce not only the tilt, but also the scat- 
ter associated with the plane. In the g-band we calculate: 
logr c = (2.12 ±0.02) log a + (0.19 ± 0.01) fj* - (7.92 ±0.07); 
and similar results in the i and z-bands. This reveals that 
the slope of the FP is i ndependent of waveleng th, analogous 
to the results found by IBernardi et all |2003ch . 

We plot the fundamental plane at different redshifts in 
Fig. 1101 The radius, velocity dispersion, surface brightness 
and mass-to-light ratios of the model galaxies all evolve with 
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Figure 9. The fundamental plane for GALFORM early-type galaxies 
(gray shading) compared with the observational data from Ber05 
(hatched shading). The shading denotes the 10 to 90 percentile 
interval. 



time (see §5), so one might expect to see some evolution in 
the fundamental plane itself, unless the changes in these 
quantities occur in such a way as to cancel out any evolu- 
tion in the locus of the plane. Ber05 report evolution in the 
position of the fundamental plane which corresponds to a 
change in the mean galaxy surface brightness of Afi e m —2z. 
We find no clear evidence for evolution in the model predic- 
tions over the same redshift interval. In section [6] we show 
the predictions for the scaling relations over a wider baseline 
in redshift. 



5 THE DEPENDENCE OF THE STRUCTURAL 
PROPERTIES OF ELLIPTICALS ON THE 
PHYSICAL INGREDIENTS OF THE MODEL 

Our calculation of the sizes of galactic spheroids contains 
several steps and is sensitive to some of the physical ingre- 
dients of the galaxy formation model more than others. The 
beauty of semi-analytical modelling is that we can switch off 
or vary particular assumptions or processes to isolate their 
impact on the model predictions. Such a study is only pos- 
sible to a very limited extent in fully numerical simulations 
of galaxy formation. Moreover, the high speed of the semi- 
analytical calculations compared with a numerical simula- 
tion allows us to examine many different variants in a short 
time. In this section, we seek to establish the sensitivity of 
our model predictions for the structural and photometric 
properties of spheroids to the composition of the model. For 
this purpose, we study the model predictions at z — and 
consider bulge dominated galaxies, i.e. those with a bulge- 
to-total luminosity ratio in the r-band of B/T > 0.8. The 
results of this section are presented in Figs. llllll2lll3l and ll4l 
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Figure 10. The evolution of the fundamental plane with redshift. 
The upper panel shows model predictions for redshifts z = 0.06, 
2 = 0.13 and z = 0.20; these are the median redshifts of the ob- 
servational samples plotted in the lower panel. Again, the shaded 
regions show the 10 to 90 percentile ranges of the distributions. 



which look, respectively, at how deviations from the funda- 
mental plane correlate with various galaxy properties, the 
Faber- Jackson relation between velocity dispersion and lu- 
minosity, the radius-luminosity relation and the fundamen- 
tal plane. 



5.1 The deviation from the fundamental plane 

As we noted in the previous section, there is some contro- 
versy in the literature regarding the source of the dispersion 
around the fundamental plane, with some authors arguing 
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that the scatter could be due to a number of causes, such 
as variations in the formation times of galaxies, metallic- 
ity trends in stellar populations, or differences in the dark 
matter content of galaxies. 

Fig. [11] shows how the deviation from the fundamen- 
tal plane correlates with various galaxy properties. Here, 
the quantity A represents the offset from the z — fun- 
damental plane predicted in the r-band, after applying 
the selection criteria to match the Ber05 SDSS sample: 
A = logn - (1.94 log a + 0.19^ e - 7.54). Fig. IliTa) shows 
that the deviation is correlated with r-band absolute mag- 
nitude for galaxies fainter than L*. This result shows that 
magnitude-limiting a sample might bias the determination 
of the fundamental plane. Surprisingly, luminous early-type 
galaxies (Af r — 51og/i < —20) exhibit no correlation with 
deviation from the funda mental plane, which is in agree- 
ment with the results of iBernardi et all (|2003bl) . For the 
remaining panels in Fig. 111! we only select galaxies with 
M r — 5 log h < —19.5 (i.e. brighter than one magnitude faint- 
wards of M*), in order to make our results comparable to 
observations. Fig. Illf bl reveals a strong anticorrelation be- 
tween the deviation and the r-band luminosity-weighted age 
of the galaxy, in the sense that galaxies which lie above 
the fundamental plane are younger. A linear fit to the dis- 
tribution reveals A = -(0.11 ± 0.03) Age + (0.64 ± 0.09). 
This strong correlation indicates that the age of the stel- 
lar population plays an important role in determining the 
position of the ga l axy in the fundam ental plane space (see 
iForbes et aUll998l ; |Pahre et al.lll999T ). We find that the fun- 
damental plane offset is als o anticorrelated with g -r colour 
(Fig. CUc)). As noted by IBernardi et"al] l|2003d ). this is 
due to the correlation between colour and velocity disper- 
sion (see also Bernardi et al. 2005). Interestingly, we see 
in Fig. Illf dl that the total stellar mass is anticorrelated 
with the deviation from the fundamental plane, for galaxies 
brighter that M r — 5 log h < —19.5. When all the early-type 
galaxies are included, then the distribution reveals a differ- 
ent picture: similar to the trend seen in panel (a), we find 
that faint galaxies show a stronger deviation from the funda- 
mental plane. Fig. Illf e) shows that the absolute metallicity 
of the stellar population is anticorrelated with a deviation 
from the fundamental plane: metal-rich galaxies are to be 
found predominately below the mean fundamental plane re- 
lation. The relation between the pseudo-specific angular mo- 
mentum, r c er, of the bulge and the FP offset is plotted in 
Fig. Illf f). revealing a weak anticorrelation between these 
quantities. There is little correlation between the deviation 
from the fundamental plane and the fraction of the total 
stellar mass formed in the last burst of star formation trig- 
gered by a galaxy merger fFig. lllf gD. which shows that the 
presence of gas in galaxy mergers does not change signifi- 
cantly the fundamental pla ne relation. This s eems to con- 
tradict the recent results of iRobertson et all (I2006T ). How- 
ever, if we also consider faint galaxies, which tend to have 
mergers containing a larger fraction of gas, then we find 
an anticorrelation between M 4 b ™' /M* ot and the deviation 
from the FP, along the lines of that seen by Robertson et al. 
The relation between the deviation from the FP and halo 
mass is shown in Fig. Hlf tO . We find an anticorrelation be- 
tween these two quantities, such that galaxies which lie in 
more massive haloes are found below the main FP relation, 
i.e. cluster galaxies should lie below the mean fundamental 



plane. In Fig. Illf i). we show that the distribution of the 
FP offset for central galaxies resembles that predicted for 
satellite early-types. 



5.2 The physics of the model and the scaling 
relations 

In this section, we examine how the predictions of the Baugh 
et al. model change if one ingredient at a time is varied. 
These variant models are not necessarily acceptable galaxy 
formation models, because they may not give as good a 
match to the local data used to calibrate the model pa- 
rameters as was the case for the fiducial Baugh et al. model. 
We also show the predictions of the Bower et al. model, as a 
further example of a variant model. In this case, many ingre- 
dients have been changed from the ones used in the Baugh 
et al. model, as explained at length in Section 2. 

The first ingredient we test is the adiabatic contraction 
model used to take into account the gravitational pull of 
the baryons on the dark matter. To recap, the condensa- 
tion of baryons at the centre of the dark matter halo pro- 
vides an additional gravitational force on the dark matter 
which causes it to move inwards, thereby increasing the den- 
sity of dark matter in the central part of the halo. This in 
turn alters the gravitational force on the baryons due to 
the dark matter. The degree of contraction is computed by 
exploiting the fact that, in a slowly varying potential, the 
action integral, <f> pidqi, is an adiabatic invariant for each 
particle of mass i, where Pi is the conjugate momentum of 
the coordinate qj (Barnes fc Whitelll984l : iBlumenthal et al.l 
1986; Ijesseit et al.ll2002t ) If we assume spherical symmetry 
and circular orbits, the action integral simplifies to the con- 
servation of angular momentum in spherical shells, rM(r). 
The adiabatic contraction of the dark matter leads to a 
more centrally peaked halo density. The main consequence 
of switching off the adiabatic contraction of the dark matter 
halo is that the half-mass radius of the spheroid increases 
(Fig. ll3f a')'). The radii of bright galaxies increase by a larger 
factor than those of faint galaxies, leading to a steepening 
of the radius-luminosity relation. The slope of the radius- 
luminosity is in much better agreement with the observed 
slope on omitting adiabatic contraction, although the model 
galaxies are too large overall (both spheroids and discs). 

Next we ignore the self-gravity of the baryons when 
computing the size and effective rotation speed of the disc 
and bulge. This also means that there is no adiabatic con- 
traction. The rotation curve of the galaxy in this case is 
set purely by the dar k matter, which is assumed to ha ve an 
NFW density profile ([Navarro. Frenk. fc Whitelll997l 'l. The 
consequences of this change are a flattening in the veloc- 
ity dispersion- luminosity relation fFig. H2f b)). with brighter 
galaxies displaying a lower velocity dispersion, and a uni- 
form increase in the radius of the spheroid (Fig. I13f b)). In 
combination, these changes result in a different projection of 
the fundamental plane which looks flatter in the projection 
which best fits the predictions of the fiducial model. 

Feedback, the regulation of the star formation rate due 
to the reheating and ejection of cooled gas following the in- 
jection of energy into the interstellar medium by supernova 
explosio ns, plays an impor tant role in setting the sizes of disc 
galaxies l|Cole et all 2000). The strength of SNe feedback is 
quantified in the model by the parameter f3 (for details see 
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Figure 11. The dependence of the deviation (defined as A = logr e — (1.941og a + 0.19/x e — 7.54)) from the fundamental plane in 
the model on various galaxy properties: (a) the r-band absolute magnitude; then, for M r — 51og/i < —19.5 galaxies (b) the r-band 
luminosity-weighted age, (c) the g — r colour, (d) the total stellar mass, (e) the stellar metallicity, (f) the pseudo-specific bulge angular 
momentum, (g) the ratio between the mass of the stars formed in the last burst, M^ UTBt , and the total stellar mass, Af* ot , at the present 
day, (h) the halo mass; and (i) whether or not the galaxy is a central galaxy or a satellite. The contours are indicative of the density 
of model galaxies. In the panel (i), the hatched histogram represents the deviation from the fundamental plane for the central galaxies, 
and the grey histogram shows the distribution for the satellite galaxies. 



ICole et alj|2000h : (5 = {V^ ot /V disc ) a ^ , where V hot and a hot 
are parameters and Vdisc is the rotation speed of the disc 
at the half mass radius. The mass of cold gas which is re- 
heated is given by M re heat = Pijj, where ip is the star forma- 
tion rate. In the Baugh et al. model, the values adopted for 
these parameters are: cVhot = 2 and Vhot = 300 km s -1 . We 
show the impact of reducing (by setting Uhot = 100 km s -1 ) 
and increasing (by setting «hot = 600 kms -1 ) the strength 
of supernova feedback in Figs. [12] 1131 and [T4l-(c)-(d). Cole 
et al. demonstrated that increasing the strength of super- 
nova feedback results in gas cooling to form stars in larger 
haloes, which leads to larger discs. Conversely, reducing the 
feedback allows gas to cool and form stars in smaller haloes 
resulting in smaller discs. These trends are reproduced in 
Figs. [T3T c) and (d). There is little change in velocity dis- 



persion on changing the strength of the supernova feedback. 
The shift in the zero-point of the radius-luminosity relation 
produces a change in the location of the fundamental plane 
(see Fig. [11(c)- (d)). 

In the fiducial GALF0RM model, spheroids are the end 
products of galaxy mergers. As we explained in § [2] the 
radius of the merger remnant is determined by conserving 
the binding energy of the individual galaxies involved in the 
merger and their relative orbital energy. The contribution of 
the orbital energy to the energy budget is parameterized by 
/orbit: the standard choice is to set / or bit = 1 and to include 
the full orbital energy in the calculation of the remnant size. 
In Figs. 1121 1131 and [141 (e). we show the effect of removing 
the contribution of the orbital energy from the calculation 
of the radius of the spheroid produced by mergers, i.e. we 
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Figure 12. The sensitivity of the model predictions for the Faber-Jackson relation to: (a) switching off adiabatic contraction, (b) 
omitting the self-gravity of the baryons, (c) reducing the strength of supernova feedback, (d) increasing the strength of SNe fe edback, 
(e) the omission of the orbital energy from the calculation of the size of the merger remnant and (f) using thelBower et al.l (120061 ) model 
with AGN feedback. In each panel, the grey line shows the median prediction from the reference model IIBaugh et al.l l2005). at z = 0. 
The black solid line shows the median for the variant model. The errorbars indicate the 10 to 90 percentile of the predictions. The dotted 
line in each panel shows the observed relation for SDSS early- type galaxies for reference ijBernardi et al.ll2003bl V 



set /orbit = 0. Perhaps surprisingly, this change results in 
an imperceptibly small change in the radius of the spheroid, 
except in the case of the brightest galaxies. 

Finally, we consider the model of iBower et al.l ([2006) , 
who implemented an AGN feedback scheme into GALF0RM, 
in which cooling flows are quenched in massive haloes at low 
redshift. As a result of this change to the cooling model in 
GALF0RM, IBower et al.l (|2006h were able to produce improved 
matches to the local B and K-band luminosity functions, 
the observed bimodality of colour distribution and the in- 
ferred evolution of the stellar mass function. In Figs. 1121 
1131 an d [TUffl we plot the scaling relations for IBower et alj 
(2006) model. Though the model performs quite well in re- 
producing the local fundamental plane of early-type galaxies 
and the Faber-Jackson relation, the radius-luminosity rela- 
tion for bright galaxies is substantially different from both 



the ob servations and from the predictions of the lBaugh et al.l 
i|2005l , model: luminous galaxies in the Bower et al. model 
are up a factor of three smaller in radius than in the Baugh 
et al. model. 

We also considered a variant of the Baugh et al. model 
in which the Kennicutt IMF was used in starbursts, in place 
of the top-heavy IMF. This produces scaling relations for 
early type galaxies which look very similar to those pre- 
sented for the Bower et al. model in Figs I12II13I and 1141 
with the main change being a shift in the predicted radius- 
luminosity relation. When the IMF is changed, the yield and 
recycled gas fraction are also changed accordingly, which af- 
fect the rate at which gas cools and alter the star formation 
timescale. This suggests that the primary difference in the 
predictions of the Baugh et al. and Bower et al. models is 
due to the choice of the IMF used in starbursts, in spite of 
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Figure 13. The sensitivity of the relation between radius and luminosity to: (a) switching off adiabatic contraction, (b) omitting the 
self-gravity of the baryons, (c) reducing the strength of supernova feedback, (d) increasing the str ength of SNe f eedba ck, (e) the omission 
of the orbital energy from the calculation of the size of the merger remnant and (f) using t he iBow er et al l j2006T ) model with AGN 
feedback. In each panel, the grey line shows the median prediction from the reference model llBaueh et alj|2005l) . at z = 0. The black 
solid line shows the median for the variant model. The errorbars indicate the 10 to 90 percentile of the predictions. The dotted line in 
each panel shows the observed relation for SDSS early-type galaxies for reference l|Bernardi et al.ll2003tJ) . 



the other differences between the models outlined in Section 
2. 



The results in this section suggest that the scaling re- 
lations of early-type galaxies are essentially insensitive to 
variations in some of the model parameters. However, this 
should not be interpreted as implying that these observa- 
tions are of limited value in constraining the models. It 
should be remembered that our starting point is a model 
of galaxy formation which has already successfully passed a 
range of comparisons with observed galaxy properties. 



6 THE EVOLUTION OF SCALING 
RELATIONS 

We now present the GALFORM predictions for the evolution 
of the structural and photometric properties of early-type 
galaxies with redshift. In this section, we consider the evo- 
lution over a much wider baseline in redshift than we ad- 
dressed in the previous section. Furthermore, in order to get 
a clear picture of the nature of the evolution, we relax some 
of the selection criteria which we applied to the model out- 
put in previous sections, where the goal was to mimic the 
Ber05 sample selection as closely as possible. The only selec- 
tion we apply in this section is that the bulge must account 
for at least 80% of the total luminosity in the rest-frame 
B-band. 

It is important to be able to disentangle changes in the 
typical stellar populations of early-type galaxies with red- 
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Figure 14. The sensitivity of the fundamental plane to: (a) switching off adiabatic contraction, (b) omitting the self-gravity of the 
baryons, (c) reducing the strength of supernova feedback, (d) increasing the st rength of SNe feed back, (e) the omission of the orbital 
energy from the calculation of the size of the merger remnant and (f) usin g thelBower et al.| d2006f) model with AGN feedback. In each 
panel, the grey line shows the median prediction from the reference model IIBaugh et al.l l2003!), at z = 0. The black solid line shows the 
median for the variant model. The errorbars indicate the 10 to 90 percentile of th e predictions. The dotted line in each panel shows the 
observed relation for SDSS early-type galaxies for reference ijBernardi et al.ll2005l V 



shift from evolution in their structural properties. Hence, 
we first examine the predicted evolution in the mass-to- 
light ratio of early-type galaxies in Fig. 1151 The stellar pop- 
ulations of early-type galaxies at z = 1 in GALF0RM have 
lower mass-to-light ratios by a factor of ~ 3 compared with 
the early-types at z = 0. This result is in agreement with 
th e change in mass-to- light ratio inferred from observations 
by Ivan de Yen et all (I2003T), d(\og M / L) / dz = -0.47 ± 0.11 
fsee also Ivan Dokkum et all l200fj ; Ivan de Wei et al.l 120051 ; 
iTreu et al.ll2005T l. We find a weak dependence of mass-to- 
light ratio on total stellar mass. However this is somewhat 
lower than what is seen obs ervationally (see IWuvts et al.l 
|2004| ; Ivan de Wei et al.ll20o"5h . Note that when we restrict 
our attention to bright galaxies (i.e. those with Mb — 
5 log ft < —19.5), the prediction for the median mass-to- 
light ratio steepens considerably, bringing the model pre- 



dictions into much better agreement with the observational 
estimates. 

The evolution with redshift of the Faber- Jackson rela- 
tion is shown in Fig. 1161 We show the correlation between 
velocity dispersion and the B-band magnitude in the rest- 
frame for the local universe (z — 0), z — 0.5 and z = 1. 
Fig. 1161 shows that the model predicts differential evolution 
in velocity dispersion with rest-frame luminosity; at brighter 
luminosities, the velocity dispersion drops by up to a factor 
of ~ 3 between z = and z = 1, whereas for fainter lumi- 
nosities, the change in velocity dispersion is much more mod- 
est. These result s are similar to the ones found in observa- 
tiona l studies (cf. ldi Serego Alighieri et al.| [2005; Red a et "all 
2006). We note that the scatter around the FJ relation seems 
to increase slightly with redshift. Furthermore, it is clear 
that the slope of the relation for faint early-type galaxies 
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Figure 15. The predicted evolution with rcdshift of the mass-to- 
light ratio in the rest-frame S-band, plotted against stellar mass. 
The grey line shows the prediction for z = 0, the short dashed line 
for z = 0.5, the solid black line for z = 1. The long dashed line 
shows the prediction for the median mass-to-light ratio at z = 1, 
when only considering galaxies with Mb — 51ogfe < —19 . 5. The 
dotted lines show the relations found by [j0recnscn et al.l (|2006h 
for the Coma cluster (grey), z ss 0, and a high redshift sample 
(black), z as 1. 



Figure 16. The predicted evolution with redshift of the relation 
between velocity dispersion and rest-frame £?-band luminosity. 
The grey solid line shows the Faber-Jackson relation in the local 
universe (z = 0), while the short-dashed line shows it at z = 0.5 
and the black solid line at z = 1. The errorbars show the 10 to 90 
percentile range of the distr ibution. The gray triangles show local 
data from lReda et al.l (2006) and the stars repres ent the sample of 
early- type galaxies in the K20 survey bv ldi Serego Aliehieri et al.l 
(2005), with 2«1. 



is shallower than that of the bright-end, which resembles 
the results for faint galaxies found by some authors (e.g. 
iMatkovic eral1l200tjl : lDavies et al.lll983l ). 

Fig. [T7] shows how the relation between radius and lu- 
minosity varies with redshift. As we saw in Section 4, the 
model predicts that the brightest early-types are too small. 
Nevertheless, it is interesting to look at the predictions at 
different redshifts. The primary agent behind the shift in 
the predictions is the passive evolution of the stellar popu- 
lations in the elliptical galaxies. Also, at a fixed stellar mass, 
the galaxy radii decrease with redshift (Coenda et al 2007; 
in preparation). As redshift increases, we see the bulk of 
the stars in ellipticals when they were younger and hence 
brighter. There is no significant trend in the size of the scat- 
ter in this relation with re dshift. 

iTruiillo et ail ((2004b) estimated the evolution of the 
radius-lu minosity relation of b right galaxies up to z ~ 3 
(see also iBouwens et al.l |2004| ) . These authors found that 
early- type galaxies, as defined by a high value of Sersic in- 
dex, show a size evolution proportional to (1 + 2: )- 101 ± 08 . 
this is comparable to the amount of evolution we predict in 
Fig. [17] for all ellipticals. 



The evolution of the fundamental plane has long been 
used to study changes in the stel l ar populations of galaxies 
(e.g. Ivan Dokkum fc Franxl Il996l; Ivan Dokkum et al.|[200ll ; 



iGebhardt et alj|2003l; 


van de Wei et al.ll2006l). As previously 


noted. iBernardi et al. 


(2003c) found evolution in the funda- 



mental plane which is consistent with the passive aging of 
the stellar population, A^i e « — 2z, but without any notice- 



able difference in the slope. There is a general consensus 
in the literature regard ing the nature of the evolution of 
the fundamental plane dGebhardt et al.ll2003l ; IZiegler et alj 
120051 : 1 J0rgensen et al.ll2006l l. However, as we have shown in 
the previous section, deviations from the fundamental plane 
relation are linked to several galaxy properties. 

In Fig. fT8 l we plot the model predictions in the 
Ijgirgensen et alj (|2006l ) projection of the fundamental plane, 
logr e = 1.2 log a + 0.33/i e — 9.1, for galaxies with Mb — 
5 log ft < —19.5. In this projection, we find no evolution in 
the slope or offset of the fundamental plane up to z = 1, for 
galaxies with log r c > 0.3, which is contrary to the claims 
made from observations. As we noted in Fig. Illf a) , the evo- 
lution found in observational studies might be partly due 
to the correlation between the magnitude and the devia- 
tion from fundamental plane, i.e. magnitude-limiting sam- 
ples might induce the zero-point of the FP to shift. On the 
other hand, the effective radius of bright early-type galax- 
ies in GALF0RM is smaller than observed (see Fig. [5]), and 
evolves with redshift, which will complicate any inferences 
drawn from the evolution of the FP. Interestingly, in the 
small radius regime, the evolution of our predicted funda- 
mental plane shows an offset similar to that observed. 



7 DISCUSSION AND CONCLUSIONS 

We p resented tests of the model proposed by I Cole et all 
(2000) to calculate the scale sizes of the disc and bulge com- 
ponents of galaxies. This is currently the most sophisticated 
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Figure 17. The predicted evolution with redshift of the rela- 
tion between radius and luminosity. The grey solid line shows the 
relation in the local universe, the black dashed line at z = 0.5 
and the black solid line at z = 1. The errorbars show the 10 to 
90 percentile ra nge of the model p redictions. The triangles show 
local data from[Rcda et al. (2006) a nd the stars are galaxies in 
the K20 survey byldT Scrego Alighicri et al ] j2005l) . with 2 Ri 1. 

model in use in semi-analytical codes to compute the radii 
of galaxies. In brief, the model assumes that galactic discs 
have an exponential profile and that spheroids follow an r 1 ^ 4 
law in projection. The hot gas atmosphere in dark matter 
haloes is assumed to have the same specific angular momen- 
tum as the dark matter. Gas is assumed to retain its angular 
momentum as it cools to form a galactic disc. The size of 
a merger remnant is computed by conserving the sum of 
the binding and orbital energies of the merging galaxies and 
applying the virial theorem. The self-gravity of the baryons 
and their impact on the distribution of dark matter in the 
central parts of the halo are taken into account. Cole et al. 
demonstrated that this model predicts scale length distri- 
butions for galactic discs which are in excellent agreement 
with observations. 

In this paper, we have carried out the first tests of 
the model predictions for the structural properties of early- 
type galaxies and the evolution of these relations with red- 
shift, using the pub lished models of lBaugh et al.l (|2005l ) and 
iBower et al.l (|2006l ) . The Baugh et al. and Bower et al. mod- 
els differ in a number of ways, as set out in Section 2. Two 
of the main differences are the manner in which the mod- 
els prevent the overproduction of bright galaxies and in the 
IMF assumed in starbursts. Bower et al. use AGN heating to 
switch off the cooling-flow in haloes with a quasistatic hot 
gas atmosphere, whereas Baugh et al. invoke a superwind 
which ejects gas that has already cooled. Perhaps contro- 
versially, Baugh et al. adopt a flat IMF in starbursts, and 
a standard solar neighbourhood IMF for quiescent star for- 
mation; in Bower et al., a standard IMF is assumed in all 
modes of star formation. We emphasize that, for the major- 



Figure 18. The predicted evolution of the fundamental plane 
with redshift for galaxies brighter than M-q — 51og/i < —19.5. 
The grey solid line shows the relation in the local universe, the 
dashed black line at z = 0.5 and the black solid line at z = 1. The 
errorbars show the 10 to 90 percentile range of t he model predic- 
tions. The dotted line shows the relation found bv |j0rgensen et al. 
[|l996h . The gray triangles represent z = data fr om Rcda et al. 
(2006) and the stars show data from Idi Serego Alighieri et al. 
120051 ). at * Ri 1. 



ity of the results presented, we have not adjusted any of the 
model parameters in order to improve the predictions for the 
fundamental plane and its projections. The one exception 
is where we exploit the modular nature of semi-analytical 
models to vary or switch off various physical ingredients of 
the model in order to assess their influence on the model 
predictions (Section 5.2). 

The model enjoys some notable successes. We demon- 
strated that the model can match the abundance of early- 
type galaxies in the SDSS sample of Ber05. We also obtain a 
reasonable match to the Faber- Jackson relation between ve- 
locity dispersion and luminosity and its evolution with red- 
shift, albeit with a shallower slope than measured by Ber05. 
Furthermore, we find a relation between velocity dispersion 
and age which is in excellent agreement with recent obser- 
vations. Perhaps most impressively, the fundamental plane 
predicted by the model is in g ood agreement with th at in- 
ferred for SDSS early-types bv lBernardi et all (|2003ch . The 
deviation from the FP relation reveals a strong correlation 
with luminosity, age, colour, stellar mass and metallicity: 
galaxies that lie above the mean fundamental plane relation 
are more luminous, younger, bluer, less massive and metal- 
poor. Furthermore, the feedback processes and clustering, 
as given by the pseudo-specific angular momentum of the 
bulge and the halo mass respectively, seem to play a role in 
defining the plane. 

Nevertheless, despite these achievements, there are 
some model predictions which disagree with the observa- 
tions. Formally, the slope of the predicted Faber- Jackson re- 
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lation is at odds with that measured by Ber05, although the 
overlap between the model galaxies and observations in this 
projection remains impressive. However, perhaps the most 
striking discrepancy is the slope of the radius-luminosity re- 
lation; the model predicts a significantly flatter radius lu- 
minosity relation than is observed. Whereas the model pre- 
dictions for the effective radii of faint spheroids are in good 
agreement with the data, the brightest galaxies are up to 
a factor of three smaller in the model. Our results suggest 
that, in the model, the brightest spheroids have less specific 
pseudo-angular momentum (i.e. jb = r^a, this is a definition 
of convenience; see Section[2]2]) than is the case for observed 
galaxies. This could be due to the model underpredicting the 
galaxy mass for a given luminosity. Somewhat surprisingly, 
the predicted slope of the radius-luminosity relation is in 
much better agreement with the observations if the adiabatic 
contraction of the halo is switched off (although, in this case, 
the model galaxies are uniformly too large without adjust- 
ing other parameters) . The adiabatic contraction of the halo 
in response to the presence of condense d baryons has been 
tested against numerical simulations (e.g.ljesseit et al.ll2002l 
ISellwood fc McGaughll2005l ; IChoi et al.ll2006t ). Our prescrip- 
tion for computing the size of merger remnants could become 
inaccurate if there is a significant fraction of mass in the form 
of cold gas. 

The other significant discrepancy is the evolution with 
redshift of the zero-point of the fundamental plane. The 
model predicts no evolution in the zero-point of the funda- 
mental plane. This is at odds with the evolution inferred ob- 
servationally, which is consistent with the shift in the mass- 
to-light ratio expected for a passively evolving stellar popu- 
lation. This discrepancy is intriguing, as the model does pre- 
dict a decline in the mass-to-light ratio of early-types with 
increasing redshift of the magnitude expected for passive 
evolution. The lack of evolution in the predicted fundamen- 
tal plane therefore points to a compensating change in one 
of the other projections; the effective radii of galaxies also 
evolve with redshift in the model. This serves as a caution- 
ary note to observational studies which interpret a shift in 
the fundamental plane in terms of a corresponding change in 
the typical mass to light ratio. The correlation between lumi- 
nosity and the deviation from the fundamental plane shows 
that part of the evolution found by observational studies 
may in fact be due to the construction of magnitude limited 
samples. 

In summary, the prescription outlined by Cole et al. 
for computing the radii of discs and bulges enjoys many 
successes, but displays a few important disagreements with 
observations. The solution of these remaining problems will 
require enhancement of the model to compute galaxy sizes, 
guided by the results of numerical simulatio ns of the growth 
of disc galaxies and galaxy mergers (e.g. lOkamoto et al.l 
120051 ; iRobertson et alj|2006r i . 
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